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ABSTRACT 

Radio astrometry of nearby, low-mass stars has the potential to be a powerful 
tool for the discovery and characterization of planetary companions. We present 
a Very Large Array survey of 172 active M dwarfs at distances of less than 10 pc. 
Twenty nine stars were detected with flux densities greater than 100 yuJy. We 
observed 7 of these stars with the Very Long Baseline Array at milliarcsecond 
resolution in three separate epochs. With a detection threshold of 500 /iJy in 
images of sensitivity Icr ~ 100 /iJy, we detected three stars three times (GJ 65B, 
GJ896A, GJ 4247), one star twice (GJ 285), and one star once (GJ 803). Two 
stars were undetected (GJ 412B and GJ 1224). For the four stars detected in 
multiple epochs, residuals from the optically-determined apparent motions have 
an rms deviation of ~ 0.2 milliarcseconds, consistent with statistical noise limits. 
Combined with previous optical astrometry, these residuals provide acceleration 
upper limits that allow us to exclude planetary companions more massive than 
3 — 6 Mjup at a distance of ~ 1 AU with a 99% confidence level. 

Subject headings: astrometry,stars:activity,stars:early-type,stars:planetary sys- 
tems,radio continuum: stars 



1. Introduction 

The study of extrasolar planets provides an important link between the study of star 
formation, the study of our own solar system, and the search for extraterrestrial life. Dur- 
ing the past two decades, radial velocity surveys have over 300 extrasolar planets with 
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13Mj (iButler et al.l l2006bl . www.exoplanets.org; www.exoplanet.eu). Doppler 
surveys show that ^10.5 % of nearby solary-type (FGK) stars have planets with orbital pe- 
riods less than 2000 days ( Gumming et al. 20081 ). 



Unfortunately, several factors make it difficult for radial velocity searches to detect 
planets around M dwarfs (e.g., paucity o f narro w spectral lines, reduced flux, photospheric 
and chromospheric activity; see I Joergend (120061)'). From the eight M star planetary systems 



detected by radial velocity techniques (IBailey et al.ll2009l). GJ 876 is p articularly interesting 
because it has a short-period ~ 7.5 planet (IRivera et al.l l2005l ) . as well as a pair of 
massive planets with o rbital periods of ^ 30 and ~ 60 days in a 2:1 mean motion resonance 



( iLaughlin et al.ll2005l ). [Benedict et al.l (120021 ) measured an astrometric perturbation by GJ 
876d of 0.25 ± 0.06mas. GJ 849 is note worthy because it ho st s a giant plane t (~ 0.8Mj) 



with an orbital period of ~ 5.1 years (IButler et al.l l2006al ). jjohnson et al.l (|2007l ) found 
evidence for at least one and perhaps two giant planets orbiting GJ 317 at separations of 
~ AU or more. Planets around GJ 849b ar id GJ 317b are estirnated to induce astrometric 
perturbations of ~ 0.6mas and ~ 0.5mas (IButler et al.ll2006al : iJohnson et al. 



2007). Most 



recently, GJ 832 was found to have a Jupiter mass planet in a nine- year orbit (Bailey et al. 



20091 ). whereas GJ 176 has a super-Earth on a nine-day orbit (iForveille et al.ll2009l ). Overall, 
the Jupiter-mass M dwarf planets have astrometric signals that are comparable to those 
detectable with radio astrometric techniques. 

The radial velocity results allow estim ates of the frequency of short-period giant plan- 
ets around M dwarfs. Endl et al. J2006h estimate a frequency of close-in Jovian planets 



around M dwarfs of < 1.3%. iJohnson et al.l (120071 ) estimate the low-mass K an dM stars 
have a 1.8 ± 1.0% planet occurrence rate, significantly less than that for solar-mass stars 
(4.2 ± 0.7%). Several groups are working to extend the Doppler technique to the near- 



infrared (e.g., iLloyd et al.l 120091 : iRamsey et al.l 120081 : iGe et al.l 120061 ). though initial results 



will still be biased towards short orbital periods. Thus, it is interesting to consider alternative 
techniques which are best suited for searching M dwarfs for planets at larger separations. 

The microlensin g technique has di scovered several extrasolar planets around low-mass 
stars. For example, ( Gaudi et ahlboosl ) present evidence for two giant planets with orbital 



separations of ~ 2.3 and ~ 4.6AU. This and two other systems (jUdalski et al.l 12005 



Dong et al.ll2008l ) provide evidence for roughly Jupiter mass planets. The host stars of the 
systems discovered by microlensing are too distant to lead to a detectable astrometric signa- 
ture. However, they do provide further evidence for a population of planets with separations 
and orbital periods that are well-suited for astrometric detection when considering nearby 
M dwarfs. Other detections jBennett et all bood : iBeauheu et aP bood : lOould et al.l l2006h 
suggest that lower-mass planets may be quite common around M dwarfs. 
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Direct ima ging has been used to detect planet candidates in wide orbits around several 
low-mass stars ( Chauvin et al .112004 . 2005 : Neuhauser et al. 2005 ). While these planets would 
induce a large astrometric perturbation, their long orbital periods would make it impractical 
to discover such planets from astrometry alone. The low density of bright M dwarfs on the sky 
makes transit searches particularly difficult. However, the potential for detecting transiting 
planets near the ha bitable zone of an M dwarf has motivated at least one dedicated M 
dwarf transit search (llrwin et al.ll2009l ). Due to the strong bias towards finding short-period 
planets, the transit technique (like the Doppler technique) is complementary to microlensing, 
direct imaging, and astrometric searches that are best suited for planets with wider orbits. 



Optica l astrometry has been used to search for planets around nearb y M dwarfs (IPravdo fc Shaklan 



19961 . I2OO3I ). leading to discovery of several low mass companions (jPravdo et al.l |2005| ). 
These studies are fore -runners of space-based astrometric planet searches such as with SIM 
Shao k Nem"^ J2009h . 



In this paper, we focus on the potential for a ground-based astrometric planet search 
using techniques of radio astronomy. Ra dio astrornetry has long been the gold standard 
for definition of celestial reference frames (IFey et al.ll2004l ) and has been used to obtain the 
most accurate geometric measurements of any astronomical technique. Astrometric results 
include measurement of the parallax and proper motion of pulsars at distances greater than 
1 kpc (IBrisken et al.l 120021). an uppe r limit to the proper motion of Sagittarius A* of a 
few km s~^ (Reid &: Brunthalerll2004l ). a < 1% distance to the Taurus star- forming cluster 



Jloinard et all boO?! ) and to Oph iuchus (Loinard et al. 2008), a ~ 5% d istance to the 



Orion Nebula star-forming cluster (jSandstrom et al.ll2007l : iMenten et al.ll2007f). and ac curate 
parallaxes to star-forming regions at distances as large as 5 kpc (e.g.. lXu et al.ll2009l ). 



Many radio-astrometric studies have been performed on stars. Very long baseline in- 
terferometry (VLBI) has been used to astrometrically discover the low-mass (M 0.1 
M(T)) companion of the southern hemisphere K dwarf AB Doradus and characterize its orbit 
(iGould et al.ll2006l ) Similarly, VLBI studies of radio-emitting stars have been used to link the 
reference frame of the HIPPARCOS satellit e to the radio extrag alactic reference frame and 
study the orbit of the ternary svster a Algol (|Lestrade et al.lll999[) . and study the structure of 



the radio emission in T Tauri stars (jPhillips et al. 



199*6^ . In fact, T Tauri itself has been the 



target of a number of astrometric studies in the radio (e.g., iLoinard et al.l 120071 ). The Very 
Long Baseline Array (VLBA) can routinely achieve an astrometric accurary of ~ 100 /zas in 
a single epoch, and it is capa ble of accuracies as high as 8 /zas under favorable circumstances 
( iFomalont fc KopeikinI 120031 ). To use this technique the target source must have a suffi- 
ciently high brightness temperature to be detected by a high resolution radio interferometer. 
For the VLBA, brightness temperatures must be > 10'' K, requiring nonthermal emission. 
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Thus, the systems which can be studied are hmited to the most active stars. 



Nonthermal stellar radio en iission has be en detected fr om many stellar ty pes (iGiidel 



20021 ). including bro wn dwarfs (IBergerl 120061) . proto-stars fiBower et al.l l2003l ). and mas- 
sive stars with w inds (Dougherty et al. 20051 ) . Radio emission from the late-type stars was 
first detected by iGary fc LinskyI (jl98ll ) and it originates in cyclotron emission due to non- 
relativistic electrons in the coronal plasma. Only late- type stars are sufficiently bright, nu- 
merous, and low mass to provide a large sample of stars suitable for large-scale astrometric 
exoplanet searches. Radio astrometric searches can determine whether or not M dwarfs, the 
largest stellar constituent of the Galaxy, are surrounded by planetary systems as frequently 
as FGK stars and how the planet mass-period-eccentricity relation varies with stellar type. 
While Doppler and transit methods can constrain this distributiuon at short-orbital periods, 
astrometry is best suited for studying planets at a few AU, where they induce a large astro- 
metric signal and it is still practical to observe the system for multiple orbital periods. The 
population of gas giants at a few AU ar ound low mass stars i s an important d iscriminant 



between planet form ation models (e.g., iLaughlin et al.l l2004l : llda fc LinI l2005l : iBosd 12006 
Kennedv et aPboO?! ). 



Radio astrometric searches for planets have a number of unique qualities. First and 
foremost, searches in the radio observe stars that because of their activity and variability 
are not good targets for radial velocity or transit studies. Thus they are highly complemen- 
tary to those carried out using other techniques. Furthemore, astrometric studies of reflex 
motion have the ability to fully characterize the orbits and masses of the detected planets, 
without the degeneracies inherent to radial velocity techniques. Finally, they are sensitive 
to long-period planets with sub- Jovian masses provided there is a long-enough time base- 
line of observations, and they naturally provide absolute astrometric positions tied to the 
extragalactic reference frame. 

The most serious limitation to astrometric accuracy may be from stellar activity that 
results in an astrophysical "jitter" added to the true source position. Most evidence, however, 
indicates tha t this jitter i s smal l enough to permit exoplanet searches around nearby stars. 
For instance, IWhite et al.l (Il994j ) model the radio emission of dMe stars as originating within 
~ 1 stellar radius of the photosphere. At a distance of 10 pc for a M5 dwarf, a stellar radius 
is ~ 0.1 mas, an order of magnitude smaller than the astrometric signature of a Jupiter 
analog. The few dwarf stars detected with VLBI appear to be compact, supporting this 
result jBenz et al.lll998l : iGuirado et al.ll2006h . 



High quality astrome tric positions of radio s t ars are also crit i cal fo r connecting radio an d 
optical reference frames (jPerryman et al. 1997 : Lestrade et al. 1999; Boboltz et al. 2007 ). 
Additionally, these observations will produce sizes, morphologies and brightness tempera- 
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tures critical for the study of p hysical processes in active stars, which are poorly understood 



[e.g., 



Gudelll2002l : lBergei]l2006f ). 



In this paper, we provide results from a flux density survey of nearby, low mass stars 
(§2 and 3). In §4, we present multi-epoch, high-resolution astrometric observations of a 
subset of the detected stars. In §5, we calculate limits on companions based on the radio 
measurements and archival optical astrometry. We summarize in §6. 

The observations described in this paper constitute a preliminary survey for the Radio 
Interferometric Planet (RIPL) search. RIPL is a program with the VLBA and the 100-m 
Green Ban k Telescope with th e goal of astrometric detection of companions to nearby, low 
mass stars (IBower et al.l 120071 ). A sample of 30 stars taken from the VLA surveys described 
below will be observed 12 times over ~ 3 years with an astrometric accuracy of ~ 0.1 mas, 
sufficient to detect Jupiter mass companions at a radius of 1 AU. The results of this paper 
demonstrate that jitter in stellar position is not a limiting factor for RIPL. 



2. Sample Definition 



St ars were draw n from two samples: active stars in the California & Carnegie planet 
search Jwrightlboosll^ and X-ray selected n earby stars which are too active for radial velocity 
searchs (NEXXUS ISchmitt fc Liefkell2004j ). X-ray fluxes are from the ROSAT all sky survey. 
The most active of the California & Carnegie stars had a lower X-ray luminosity than the 
least active of the NEXXUS stars. All stars are M d warfs at D < 10 pc. In addition , 



we include three star s, GJ436 rM2.5 iButler et al.l 120041 ) . e Eri rK2V ICampbell et al.lll988 



Gumming et al.lll999l ) and HD 131156 (G8V), which are chromospherically active, nearby, 
and show ra dial velocity signatur es for planetary companions (e Eri has also a nearly face-on 
debris disk; ICreaves et aD Jl998h l We also includ e the dMe star AU Mic (GJ 803) which 
has a recently discovered nearly edge-on debris disk (IKalas et al.ll2004j ). since debris disks are 
believed to be signposts of planet formation. Basic stellar data for 172 stars are tabulated 
in Table [H These data comprise extended Gliese catalog number, coordinates, error ellipse 
in position, proper motion, error ellipse in proper motion, annual parallax and its error, 
spectral type, apparent B and V magnitudes and X-ray luminosity. 



3. Flux Density Observations 



We used the Very Large Array to observe our sample in June through September 2005. 
Observations were made at 5 GHz in standard continuum mode. The array was in B, C, and 
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D configurations, giving a resolution that ranged from ~ 1 to 10 arcsec. Eacli source was 
observed for 10 minutes per epoch; some sources were observ ed multiple tim es. Standard 
calibration and imaging techniques were applied with MPS JCreisenl I2OO3I 1 The typical 
image rms was ~ 50/iJy. 

In Table [2] we list flux densities for sources detected in at least one epoch. We made 
40 detections of 29 individual stars. In Table [3l we list upper limits to the flux density for 
sources that were not detected. 

We plot the radio and X-ray luminosities for the detected and non-detected stars in Fig- 
ure [H X-ray luminosities are from the ROSAT survey and so are non-contemporaneous with 
radio obser vations . Nev ertheless, we see rough agreement with the radio-X-ray luminosity 
correlation (lGiidelll2002l ). We also separate late and early type stars at type M5 but find no 
difference in detection rates or flux densities. 



4. VLBA Astrometric Survey 

4.1. Observations and Data Analysis 

In Spring, 2006, we studied the astrometric stability of seven stars with the VLBA. 
These stars were selected as among the brightest stars from the VLA sample. For each star, 
three VLBA epochs were spread over 11 days or less (Table Hj). Observations were obtained 
at a frequency of 8.4 GHz with a recording bandwidth of 256 Mb/sec. Phase-referenced 
observations were obtained for each of the stars with an integration time on the star of ~ 1 
hour distributed over 4 to 8 hours, achieving an rms sensitivity of ~ 100 /iJy. Beam sizes 
and rms flux densities for each epoch are in Table |H 

Data were processed using the VLBA pipeline in AIPS. Standard amplitude and phase 
calibration techniques were employed. The positions and flux densities of primary and 
secondary calibrators are listed in Table [51 In the case of primary calibrators, the positions 
are those assumed for correlation. These positions were taken from VLBA calibrator lists 
and are typically accurate in an absolute sense to a milli-arcsecond. Secondary calibrator 
positions were determined by phase referencing to the primary calibrator. The reported 
positions and flux densities are the averages over the three epochs. 
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4.2. Astrometry 

Of the 7 stars observed by the VLBA, three stars (GJ 4247, GJ 65B, and GJ 896A) 
were detected in all three epochs, one star (GJ 285) was detected in two epochs, one star 
was detected in only one epoch (GJ 803), and two stars were not detected in any epoch (GJ 
1224). The detections had flux densities ranging from 0.5 to 3.9 m Jy (equivalent to SNR~ 
5 to 20). Images were reasonably well-fit as point sources; we discuss possible deviations 
from point source images later. Best-fit positions and flux densities are given in Table [61 We 
show images of the stars and their primary and secondary calibrators in Figures E], [3l HI [5l 

El El El 

In the second epoch of observation, G J 285 was not detected and we show an image at its 
expected position. The secondary calibrators were imaged indicating that phase-referencing 
was successful. The absence of GJ 285 thus must be attributed to a flux density below the 
3cr-limit of ~ 300 /iJy or to resolved structure. 

GJ 65B has a large apparent motion (~ 16 mas/day) and was significantly variable 
during the experiment. The apparent motion due to both its large proper motion and large 
parallax causes the source to move more than a beam width during the 5 hours of observation. 
We split epochs into two time ranges and found significantly variability in the flux density. 
In the first epoch, GJ 65B is detected in both halves of the epoch. In the second and third 
epochs, GJ 65B is detected in only the second half of each epoch. Observational results are 
tabulated for these four half epochs in which GJ 65B is detected. 

We detected GJ 803 in only epoch. In the second epoch, phase stability was very poor 
such that the images of the calibrators were severely distorted. In the third epoch, the 
calibrator images were of good quality but there was no detection of GJ 803 at a level of 
5cr within the field of view anticipated given the detection in the first epoch and optical 
astrometry. Imaging quality is poor for this source due to its low declination. For second 
and third epochs, we show images at the expected locations of the star, assuming optical 
proper motions and parallaxes and the position detected in the first epoch. 

Detections were not obtained in any epoch for two stars, GJ 412B and GJ 1224. Ob- 
serving conditions for GJ 412B were good on all three epochs but no source was detected. 
Compact calibrators JlllO-l-44 and J1058-I-43 were detected with consistent fiux density, 
structure, and positions in all three epochs. GJ 412B is separated from the calibrator by 
approximately 30', which is less than the separation of the secondary calibrators. Thus, poor 
phase calibration is unlikely to be the cause for the lack of detection of GJ 412B. Based on 
optical astrometry, the apparent position of GJ 412B was changing during these epochs at a 
rate ~ 0.6 mas h~^, which means that it moved more than a synthesized beam width during 
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the course of the observations. If uncorrected, this leads to an approximately 40% reduction 
in the observed peak flux density. We corrected for this effect by making images from short 
time segments, shifting the positions of these segments to account for the optical astrometric 
model, and then co-adding these segments. The star was not detected in the final co-added 
images. We searched images that were ~ 2.2" on a side, centered on the optical position. 
The likely cause of the lack of detection is low flux density. In two out of three VLA epochs, 
GJ 412B was not detected above 200 /iJy. 

In the case of GJ 1224, phase-referencing failed. This is due to the large offset in 
declination (~ 3 degrees) between the phase calibrator J1753-1843 and the star. The image 
of secondary calibrator J1809-152, which is within 0.7 deg of GJ 1224, showed the effects of 
poor phase calibration. On the other hand, the secondary calibrator J1825-1718 is within a 
degree in declination of the primary calibrator and was observed to have a relatively compact 
structure. Self-calibrated observation of J1809-152 indicate a compact source with 60 mJy 
flux density, suitable as a phase calibrator for future observations of GJ 1224. 



Discussion 



5.1. Survey Images 

The VLBA images in our exploratory survey, with the exception of that for GJ65B, are 
predominantly compact and point-like. This is consistent with the radio emission originating 
from points close to the stellar photosphere, which has a characterist ic radius of ~ O.li?© 
(equivalent to 6^, ^ 100 /ias at a distance of 5 pc; iBeatty et al.l 120071 ). 



The image of GJ65B (also known as UV GET B) is extended, as was previously seen by 
Benz et al.l (119981). The ima g e app ears to consist of two components, possibly with a variable 
position angle. iBenz et al.l (119981 ) argue that these two components are magnetic loops in 
the corona and that the star is located betwee n these two. We find t hat the flux density of 



GJ65B is variable, as was previously observed (ILinsky &: Garylll983l : iPallavicini et al.lll985 
Jackson et al.lll987l ). 



GJ 896A (EQ Pe g) has been previously detected with VLBI and found to be com- 
pact (IBenz et al.l Il995l ). GJ 285 (YZ GMi) has been previously detected as compact an d 
marginally resolved with VLBI in different epochs (IBenz fc Alellll99ll : lPestalozzi et al.ll2000l ). 
GJ 803 and 4247 have not been previously imaged with VLBI to our knowledge. We leave 
a detailed investigation of stellar activity to another paper. 
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5.2. Astrometric Results 

For the four sources (GJ 896A, GJ 4247, GJ 65B, and GJ 285) detected in multiple 
epochs, we plot their positions as a function of time in Figures |9] through [T2l Positions are 
referred to the expected position for the first epoch from the existing optically-determined 
astrometry, proper motion, and parallax listed in Table [H The se originate from Hipparcos 



observations for GJ 285 and GJ 896A (IPerrvman et al.lll997D. from a comparison of the 



Luyten and 2MASS catalogs for GJ 65B (ISalim fc Gouldl 120031) . and f rom a comparison 



of the Luyten and Tycho-2 catalogs for GJ 4247 (ILepine fc Shard 120051 ) . In Table [6], we 
tabulate the radio coordinates relative to the first epoch VLB A position (Aa, A6) and to 
the expected optical positions (Aooptj A(5opt). 

We have compared the measured positions of the radio stars with their predicted opti- 
cal positions. The uncertainties in the radio-optical position offsets are dominated by the 
relatively large astrometric uncertainties in the optical measurements, primarily due to the 
faintness of the stars and the long time baseline between the optical and VLBA measure- 
ments. As an ensemble, the mean positions are consistent with no offset. The reduced 
is 0.8 and 0.7 for a and S, respectively. The positions measured by the VLBA are consis- 
tent with the optical astrometry for GJ 4247, GJ 803, and GJ 285 in both coordinates. In 
the case of GJ896A, there is a marginal detection (3.2(j) of an offset in the absolute right 
ascension. In the case of GJ65B, the total radio-optical offset is also marginally significant. 
As we discuss below, this source is part of a binary, which complicates the interpretation of 
these short observations. 

The RMS residuals relative to the optical astrometry in each coordinate are ~ 0.1 to 
~ 0.2 mas (Figures [IDl [HI and [T^. This RMS is consistent with statistical errors in 
the VLBA observations, indicating that there is not a significant contribution to the source 
positions from stellar activity. The absence of a stellar activity contribution to the astrometry 
is a central result of this paper. It remains possible that the emission originates from a region 
that is a few stellar radii in scale but has a stable centroid on a time-scale of days. It is also 
possible that there is longer-term stellar activity that will corrupt astrometric accuracy. 

We have also compared the relative apparent motions of our sources. For all sources, 
we compute the difference in the apparent motions determined by VLBA and that from the 
optical catalogs (A/Xq, and A/x^; Table [7]). These apparent motions include the effects of 
parallax and proper motion. We then compute the implied acceleration in each coordinate 
using the time difference between the epoch of the VLBA observations and the epoch of 
the optical position determinations. For Hipparcos observations, the epoch of the optical 
observations is 1991.25. For the other two catalogs, we take the epoch of proper motion as 
the mean of the Luyten catalog epoch (1950) and the relevant modern optical catalog, either 
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2MASS (epoch 1999) or Tycho-2 (epoch 1991.25). The uncertainty in acceleration has a 
significant component set by the relatively short time-baseline of the VLBA measurements. 
For the cases of GJ 65B and GJ 4247, the VLBA observations span only 3 days, which leads 
to an error in apparent motion of ~ 0.2 mas/3 days ~ 25 mas/year. Observations that span 
a year will be one to two orders of magnitude more sensitive to acceleration. 

For all stars but GJ 65B, there is no detection of an apparent acceleration relative to the 
optical astrometry. In the case of GJ 65B there is a significant offset {6a) in the declination 
proper motion. The right ascension proper motion is fully consistent with statistical errors 
in the VLBA measurements. As we discuss below, the apparent acceleration for GJ 65B is 
consistent with the expectations for the binary orbit. 



5.3. Limits on Companion Mass and Semi-Major Axis 

For the stars for which we find an upper limit on the acceleration, we compute limits 
on companion mass and semi-major axis, assuming circular orbits. We assume a mass of 
0.1 Mq for the star. Our calculation determines the fraction of systems with which a 3a 
detection is made for a given set of planetary mass and semi-major axis. 

The refiex motion, R, of a star of mass M due to an orbiting planet of mass Mp orbiting 
a distance r away from the center of mass of the system can be simply expressed as 



Assuming a circular orbit, the instantaneous star speed and acceleration are 



dR 
'dt 
d^R 



oc 



oc 



-/Mr 



(2) 
(3) 



where r is the orbital semimajor axis of the planet and we assume m -C M. At small semi- 
major axes, limits on companion mass and semi-major axis are set by the maximum angular 
displacement of the star during the period covered by the observations. At these small 
separations, accelerations larger than the minimum observable acceleration can occur, but 
these do not produce an angular displacement detectable by the VLBA. We have calculated 
the acceleration of the star, projecting the systems over the full range of stellar positions and 
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orbital inclination angles and sum over all systems that produce a detectable offset. This 
effect then incorporates the loss due to systems with an orbital period much less than the 
baseline time. The corresponding limit on companion masses then scales as r~^. 

At large star-planet separations (r) — and therefore long periods — the limits are set 
by the acceleration due to the companion. The fraction of systems that have detectable 
accelerations for a given radius are a function of the inclination angle of the binary and the 
position angle of the two proper motion measurements. To determine our confidence limits we 
compute the acceleration for a face-on system and then project the orbit onto a full range of 
position and inclination angles and sum over all systems that produce an acceleration greater 
than our acceleration limits. The calculation does not take into account the degenerate 
condition that occurs when the period is very close to the separation between the optical 
and radio epochs. For stars with Hipparcos data, this peak is at 15 years, corresponding to 
3 AU for the 0.1 Mq star; for the other stars, the time baseline is > 30 to 50 years. The 
mass limits scales with distance to the star as (see Eq. [3]). Acceleration limits are three 
times the quadrature sum of the errors in the two coordinates. 

For the data presented here the acceleration limits correspond to a minimum in mass at 
an approximate semi- major orbital axis rmin ~ 1 AU (Figures [T3l) . That is, these observations 
are most sensitive to planets orbiting at that distance from the primary. Near this minimum 
the orbital periods are ~ 3 years and the minima in mass (Mp^mm) are in the range 2 to 5 
Mjup (Table E]). We also give the planet mass detection thresholds at 0.3 AU (Mp o.s) and 3 
AU (Mp s); these are in the range of 10 to 20 Mj and 20 to 40 Mj, respectively. 

To be conservative, we treat the detection of acceleration for GJ65B as an upper limit 
as well. For the nominal acceleration detection, a = 0.0192 ± 0.0038AU/yr^, solutions for 
r > 1.4 AU fall on the curve 



4.5MjupJ \3AU 



) ■ (4) 



GJ65B has an M dwarf companion with an orbita l period P = 26.5 2 yr, inclination 



i = 127°, and semimajor axis a = 1.95 arcsec or 6.3 AU (IGeyer et al.lll988l ). The maximum 
apparent acceleration for this system is 0.083 AU/yr^. Observations at non-optimal epochs 
(as we have in this case) will lead to reduced apparent acceleration, consistent with the 
acceleration we have observed. 

It is difficult to model exactly the contribution of the binary orbit to the proper motion 
and acceleration of G J65B given uncertainty over exact observing epochs and methods for cal- 
culating the published proper motions based on optical observations. But it appears that the 
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binary orbit was not taken into account in the proper motion calculation of ISalim fc Gould 



( 120031 ) ■ For epoch 1950 of the Luyten catalog and epoch 2006 of the VLBA observations, 
the star was near apastron while for epoch 1999 of the 2MASS catalog the star was near 
periastron. We estimate an additional error of (10,30) mas/y in the right ascencion and 
declination proper motions based on the stellar orbit. This has the effect of reducing the 
significance of A/i^ and the acceleration as to less than Sa. We also note that observations 
at low declination are prone to more systematic error in the declination coordinate. 

GJ896A is also in a close binary system with another low mass star. We are unaware 
of an orbital solution for the system and therefore cannot place limits on acceleration due 
to the companion. 



6. Conclusions 

We measured flux densities of a sample of X-ray selected low mass stars in the stellar 
neighborhood. We detected 29 of these stars, consistent with the expectations of the radio- 
X-ray correlation. Of these stars, we observed 7 with the Very Long Baseline Array and 
detected 5. Astrometry of these stars indicates that we are not limited by jitter in the stellar 
position at a level of ~ 0.2 milliarcseconds. Provided that there is not longer-term evolution 
of the radio activity of these stars, our results indicate indicate that radio-monitoring of 
these stars can be effective for detection of Jupiter mass planets. 

A comparison of radio and optical astrometry allows us to place upper limits on com- 
panions. We exclued companions of a ~ 3 to 6 Jupiter masses at a radius of ~ 1 AU for 
four stars. At radii of 0.3 and 3 AU, limits on companion masses are in the range of 10 to 
20 Mj and 20 to 40 Mj, respectively. The short time-baseline of the VLBA measurements 
limits the acceleration accuracy. Longer time-scale analysis of these sources and others will 
place much stricter constraints on companion masses. 

Additionally, several improvements to observing and analysis methods can improve the 
quality of results. Due to their large proper motion, some of these stars change position by 
amounts comparable to or larger than the size of the VLBA synthesized beam during a few 
hours of observation. Use of data integrated throughout an entire interferometric track thus 
results in degraded image quality and astrometry. We made a first-order attempt to correct 
this by splitting some observations into two segments and analyzing them independently; 
however, more sophisticated approaches in the visibility domain will generate higher quality 
results. Accurate handling of known orbits for binary stars is also necessary to achieve full 
sensitivity to planetary companions. Fortunately, known binary companions have periods 
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that are long compared to the observing period and so residuals are likely to appear as linear 
terms that can be removed through proper motion fitting. Further, we have performed a 
simple analysis of orbital parameters given the limited nature of this data; future model- 
ing must solve for orbital parameters in a more complete and thorough manner. Finally, 
greater sensitivity through increased recording bandwidth and the use of larger apertures 
can increase SNR and push astrometry accuracy to 0.1 milliarcseconds and smaller. 

This paper is the first result from the VLBA+GBT Radio Interferomctric Planetary 
(RIPL) survey. Future papers will investigate the dctectability of other stars at high angular 
resolution, the long-term stability of stellar positions, and the detection of planets through 
astrometric means. 

We thank Jason Wright for generously sharing data that helped define the RIPL target 
sample. We wish to thank Andrew Howard and Stephen White for their useful comments. 
The National Radio Astronomy Observatory is a facihty of the National Science Founda- 
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Jet Propulsion Laboratory, California Institute of Technology, under contract with the Na- 
tional Aeronautics and Space Administration. This research has made use of the SIMBAD 
database, operated at CDS, Strasbourg, France. 
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Fig. 1. — Radio and X-ray luminosities. Results are plotted separately for late type (filled 
symbols) and early type (unfilled symbols), where the dividing spectral type is M4. Squares 
represent detections and triangles represent upper hits. The solid line is the nominal cor- 
relation Lji = lO^^^Lx Hz. Dashed lines represent the standard uncertainty of an order of 
magnitude in the scale factor. A histogram of detected fraction as a function of Lx is also 
given. 
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Fig. 2. — Images of GJ 65B and its calibrators J0142-1633 (secondary) and J0132-16 (pri- 
mary) from all observing epochs. Two images of G J65B from the first and second half of the 
first epoch are shown. Contours are -3, 3, 4.2, 6, 8.4, 12, 16.8 times the image rms noise for 
the star and the secondary calibrator. Contours are -5, 5, 7, 10, 14, 160 times the image 
rms noise for the primary calibrators. The synthesized beam is shown in the lower left. The 
epoch is written in the upper left. 
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Fig. 3. — Images of GJ 285 and its calibrators J075H-152 (secondary) and J0739+0137 
(primary). Contours and labels described for Figure [3l 
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Fig. 4.— Images of the calibrators for GJ 412B: Jl 110+4403 (secondary) and J1108+4330 
(primary). Contours and labels described for Figure [2J 




Fig. 5. — Images of GJ 803 and its calibrators J2042-3152 (secondary) and J2056-3208 
(primary). Contours and labels described for Figure El 
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Fig. 6.— Images of GJ 896 and its calibrators J2334+2010 
(primary). Contours and labels described for Figure El 



(secondary) and J2328+1929 
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Fig. 7.— Images of the calibrators for GJ 1224: J1809-1520 and J1825-1718 (secondary) 
and J 1753- 1843 (primary). Contours and labels are as described for Figure O 
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Fig. 8.— Images of GJ 4247 and its calibrators J2203+2811 (secondary) and J2205+2926 
(primary). Contours and labels described for Figure El 
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Fig. 9. — Motion of GJ 65B in RA and Dec. Plotted radio positions are relative to optical 
position at first epoch. The sohd line shows the predicted optical position due to parallax and 
proper motion. The lower panels show the residual after removal of the optical predictions 
and the mean offset between radio and optical. 
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Fig. 10.— Motion of GJ 285. See Figure [9] for details. 
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Fig. 11.— Motion of GJ 896A. See Figure 1 for details. 
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Fig. 12.— Motion of GJ 4247. See Figure [9] for details. 
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Fig. 13. — Region of planetary mass and radius phase-space rejected by acceleration upper 
limits. Contours are for parameters for which 99%, 90%, 50%, and 10% of systems would be 
detected with Sa confidence. Parameter space above the curves is rejected. Wiggles in the 
curves are due to logarithmic gridding of the model. 
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1.66, 1.22, 67 


549, 


,30 


1 


,58 


M4Ve 


11, 


,28 


9.54 


25.85 


701 


18 


05 07.58 


-3 -1 -52, 


,75 


10.60,08.53, 90 


570.10, -332, 


,50 


1.22, 1.00, 90 


128, 


,20 


1, 


,44 


Ml 


10, 


,89 


9.37 


27.13 


1224 


18 


07 32.85 


-15 -57 -47, 


,10 




-617.00, -342, 


,00 


5.00, 5.00, 176 


132, 


,60 


0, 


,00 


M4.5 


15, 


,44 


13.64 


27.93 


4053 


18 


18 57.40 


66 11 32, 


,40 


2000.00,2000.00, 45 


464.00, -469, 


,00 


100.00, 100.00, 45 


137, 


,30 


5, 


.30 


M4.5 


15, 


,29 


13.46 


27.33 


4063 


18 


34 36.65 


40 07 26. 


,44 




57.00, -203, 


,00 


2.00, 2.00, 90 


138, 


,00 


40, 


,00 


M3.5 


12, 


,84 


11.42 


26.92 


1230A 


18 


41 09.86 


24 47 13. 


,90 


2000.00,2000.00, 2 


517.00, 55, 


,00 


100.00, 100.00, 45 










M4.5 


14, 


,00 


12.40 


27.94 


1230B 


18 


41 09.86 


24 47 13. 


,90 


2000.00,2000.00, 2 


517.00, 53, 


,00 


100.00, 100.00, 45 










M5 


16, 


,10 


14.40 


27.94 



Table 1 — Continued 



GJ 


RA 


Dec 


Pos. Err 


PMa,PMi 




rr rr cw 


Sp 






Log Lx 




(J2000) 


(J2000) 


(mas,mas,deg) 


(mas/yr) 


(mas/yr, mas/yr, dcg) 


(mas) (mas) 








(erg s-1) 



725A 


18 


42 


46, 


.69 


59 37 49, 


.42 


32.44,27.69, 168 


-1327.00, 1801, 


.90 


3.61, 3.06, 164 


280, 


.20 


2, 


.57 


M3V 


10, 


.45 


8, 


.91 


26.65 


725B 


18 


42 


46, 


.90 


59 37 36, 


.65 


110.99,102.93, 162 


-1393.00, 1845, 


.50 


12.19, 11.32, 153 


284, 


.40 


5, 


.01 


M3.5 


11, 


.28 


9, 


.69 


26.65 


729 


18 


49 


49, 


.36 


-23 -50 -10. 


.44 


19.73,12.25, 94 


637.57, -192, 


.40 


2.23, 1.44, 94 


336 


.40 


1, 


.82 


M3.5 


12, 


.55 


10, 


.95 


27.78 


747A 


19 


07 


43, 


.00 


32 32 41. 


.30 




1235.00, 1130. 


.00 


5.00, 5.00, 125 


123, 


.00 


4, 


.00 


M3 


13, 


.30 


11, 


.70 


27.08 


747B 


19 


07 


45, 


.00 


32 32 54. 


.00 


18043.00,18043.00, 138 


1213.00, 1095. 


.00 


25.00, 25.00, 48 










M5 






12, 


.10 


27.08 


752A 


19 


16 


55, 


.26 


05 10 08, 


.05 


11.77,07.86, 87 


-578.80, -1331. 


.00 


1.33, 0.91, 89 


170, 


.20 


1, 


.37 


M2.5 


10 


.63 


9, 


.13 


26.94 


752B 


19 


16 


57, 


.66 


05 09 00, 


.40 


2000.00,2000.00, 179 


-592.00, -1400, 


.00 


100.00, 100.00, 45 


170, 


.20 


0, 


.00 


M8V 


19, 


.42 


17, 


.30 


27.05 


1245A 


19 


53 


54, 


.48 


44 24 53, 


.30 




443.00, -581, 


.00 


20.00, 20.00, 112 


215, 


.00 


3, 


.00 


M5.5 


15, 


.31 


13 


.41 


27.47 


1245C 


19 


53 


54, 


.90 


44 24 54, 


.00 


3000.00,3000.00, 45 


. . . , . 




. . . , . . . , ... 


220, 


.20 


0, 


.00 


M5.5 






17, 


.50 


27.47 


1245B 


19 


53 


55, 


.14 


44 24 54. 


.10 




443.00, -581, 


.00 


20.00, 20.00, 112 


206, 


.00 


3, 


.00 


M5.5 


15, 


.97 


13, 


.99 


27.47 


791.2 


20 


29 


48, 


.41 


09 41 19. 


.30 


2000.00,2000.00, 


704.00, 92, 


.00 


100.00, 100.00, 94 


112, 


.00 


3, 


.00 


M4.5V 


14, 


.69 


13, 


.04 


27.89 


793 


20 


30 


32, 


.05 


65 26 58, 


.42 


09.95,07.43, 29 


443.24, 284. 


.03 


1.15, 0.87, 29 


125, 


.60 


1, 


.11 


M2.5 


12, 


.00 


10, 


.44 


27.77 


803 


20 


45 


09 


.53 


-31 -20 -27, 


.24 


14.87,06.92, 65 


280.37, -360, 


.00 


1.68, 0.79, 65 


100 


.50 


1, 


.35 


MlVe 


10 


.05 


8, 


.61 




809 


20 


53 


19, 


.79 


62 09 15, 


.81 


07.04,05.18, 30 


1.08, -774, 


.20 


0.80, 0.62, 28 


141, 


.90 


0, 


.77 


M0.5 


10, 


.03 


8, 


.54 


27.4I7 


829 


21 


29 


36, 


.81 


17 38 35, 


.85 


24.71,09.75, 86 


1008.00, 376, 


.20 


2.91, 1.12, 85 


148, 


.20 


1, 


.85 


M3.5 


11, 


.96 


10, 


.35 


CO 


831A 


21 


31 


17, 


.80 


-9 -47 -25, 


.00 




. . . , . 




. . . , . . . , ... 










M4.5 






12, 


.60 


27./^5 


4247 


22 


01 


13 


.12 


28 18 24. 


.86 




372.11, 36, 


.48 


3.51, 2.59, 23 


111, 


.50 


3, 


.19 


M4 


13, 


.90 


11, 


.99 


28.36 


849 


22 


09 


40, 


.35 


-4 -38 -26, 


.62 


24.38,11.58, 93 


1134.90, -19, 


.71 


2.75, 1.30, 91 


113 


.90 


2, 


.10 


M3.5 


11 


.94 


10 


.42 




4274 


22 


23 


06 


.99 


-17 -36 -26, 


.40 




248. OU, -895 


.00 


231.00, 64.00, 180 


134, 


.10 


5 


.60 


M4.5 


15 


.09 


13 


.25 


27.80 


860A 


22 


27 


59, 


.47 


57 41 45, 


.15 


27.09,25.47, 64 


-870.30, -471, 


.20 


3.06, 2.91, 48 


249, 


.50 


3, 


.03 


M3 


11, 


.24 


9, 


.59 


27.72 


860B 


22 


28 


00, 


.34 


57 41 44, 


.40 


2000.00,2000.00, 85 


-711.00, -320, 


.00 


100.00, 100.00, 80 


249, 


.50 


0, 


.00 


M4 


12, 


.90 


10, 


.30 


27.72 


866 


22 


38 


33 


.62 


-15 -17 -59, 


.20 




2314.00, 2295. 


.00 


5.00, 5.00, 13 


300 


.00 


5, 


.00 


M5.5 


14, 


.14 


12, 


.18 


27.23 


867B 


22 


38 


45, 


.31 


-20 -36 -52. 


.10 




427.00, -61, 


.00 


7.00, 7.00, 68 










M3.5 


13, 


.04 


11, 


.43 


29.31 


867A 


22 


38 


45, 


.58 


-20 -37 -16. 


.08 


108000.00,108000.00, 255 


450.59, -79, 


.86 




115, 


.70 


1, 


.50 


Ml. 5 


10, 


.57 


9, 


.07 


29.31 


873 


22 


46 


49, 


.73 


44 20 02, 


.37 


11.83,09.48, 149 


-704.60, -459, 


.40 


1.40, 1.10, 148 


198, 


.00 


2, 


.05 


M3.5e 


11, 


.45 


10, 


.09 


28.99 


876A 


22 


53 


16, 


.73 


-14 -15 -49, 


.32 


32.97,13.91, 95 


960.31, -675, 


.60 


3.77, 1.58, 97 


212, 


.60 


2, 


.10 


M4 


11 


.77 


10, 


.17 


26.49 


880 


22 


56 


34, 


.81 


16 33 12, 


.36 


13.63,07.01, 31 


-1033.00, -283, 


.30 


1.60, 0.80, 30 


145, 


.20 


1, 


.22 


M1.5V 


10, 


.17 


8, 


.66 


27.09 


896A 


23 


31 


52, 


.18 


19 56 14, 


.13 


19.96,11.30, 60 


554.40, -62, 


.61 


2.25, 1.28, 56 


160, 


.00 


2, 


.81 


M3.5 


11, 


.51 


10, 


.32 


29.06 


896B 


23 


31 


52, 


.56 


19 56 13. 


.90 




602.00, 17, 


.00 


14.00, 9.00, 


160 


.00 





.00 


M4.5 


14, 


.40 


12, 


.40 


29.06 


905 


23 


41 


55, 


.01 


44 10 38. 


.90 




100.00, -1594. 


.00 


5.00, 5.00, 173 


315, 


.00 


2, 


.00 


M5.5V 


14 


.19 


12, 


.28 


27.10 


1289 


23 


43 


06, 


.28 


36 32 14. 


.00 


2000.00,2000.00, 45 


930.00, -136 


.00 


100.00, 100.00, 










M4 


14, 


.27 


12, 


.67 


27.69 


4360 


23 


45 


31, 


.27 


-16 -10 -19. 


.30 


2000.00,2000.00, 177 


-395.00, -558, 


.00 












M5 


14, 


.80 


14, 


.50 


27.66 


908 


23 


49 


12, 


.53 


02 24 04. 


.40 


19.28,06.30, 81 


995.31, -968, 


.40 


2.33, 0.77, 82 


167, 


.50 


1, 


.49 


Ml 


10 


.46 


8, 


.98 


27.12 
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Table 2. Detected Stars 



GJ Epoch Flux Density 

(YYYYMMDD) (^Jy) 



53B 


20050825 


350 ± 86 


65B 


20050807 


4271 ± 66 


84 


20050807 


196 ± 57 


102 


20050807 


182 ± 49 




20050825 


< 204 


109 


20050612 


131 ±40 




20050709 


< 138 


412B 


20050612 


152 ± 45 




20050711 


< 189 




20050903 


1300 ± 66 


557 


20050612 


236 ± 68 


644C 


20050816 


< 171 




20050824 


224 ± 57 


661AB 


20050816 


536 ± 52 




20050822 


< 147 


686 


20050609 


122 ± 39 




20050711 


< 147 


729 


20050816 


261 ± 59 




20050822 


220 ± 52 


747A 


20050809 


450 ± 133 




20050824 


< 255 


803 


20050609 


1232 ± 50 




20050612 


368 ± 53 




20050711 


222 ± 59 




20050729 


324 ± 49 


866 


20050709 


158 ± 46 


867B 


20050824 


376 ± 64 


873 


20050824 


546 ± 90 


896A 


20050807 


567 ± 52 




20050825 


1027 ± 62 



-37- 



Table 2 — Continued 



GJ 


Epoch 


Flux Density 




(YYYYMMDD) 


(/^Jy) 


1005A 


20050807 


< 204 




20050825 


348 ± 96 


1116AB 


20050903 


1482 ± 82 


1207 


20050816 


780 ± 63 




20050824 


195 ± 58 


1224 


20050816 


1271 ± 69 




20050824 


185 ± 59 




zuuuuouy 


^ 1 71 




20050824 


199 ±61 


2066 


20050709 


364 ± 55 




20050711 


< 177 


3146 


20050809 


< 267 




20050825 


281 ± 71 


3789 


20050816 


6203 ± 105 




20050824 


3139 ± 65 


4063 


20050822 


< 150 




20050816 


219 ±53 


4247 


20050816 


3054 ± 75 




20050822 


791 ± 66 


4360 


20050807 


2174 ± 130 




20050825 


< 300 
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Table 3. Undetected Stars 



GJ Epoch Flux Density Limit 

(YYYYMMDD) (/xJy) 



15AB 


20050612 


< 123 




20050709 


< 144 


48 


20050612 


< 132 




20050709 


< 147 


54.1 


20050807 


< 567 


83.1 


20050809 


< 243 




20050825 


< 210 


105C 


20050807 


< 201 




20050825 


< 231 


144.0 


20050709 


< 144 




20050711 


< 105 


166C 


20050807 


< 618 


169.1A 


20050809 


< 234 




20050825 


< 231 


176 


20050709 


< 216 




20050711 


< 255 


205 


20050709 


< 153 




20050711 


< 144 


226 


20050612 


< 141 




20050709 


< 156 


251 


20050722 


< 180 


273 


20050709 


< 156 




20050711 


< 162 


338A 


20050709 


< 156 




20050711 


< 120 


382 


20050722 


< 147 


388 


20050903 


< 567 


393 


20050722 


< 156 


406 


20050903 


< 735 


408 


20050612 


< 126 
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Table 3 — Continued 



GJ Epoch Flux Density Limit 

(YYYYMMDD) (/xJy) 





20050722 


< 141 


411 


20050612 


< 138 




20050722 


< 186 




20050903 


< 192 


412A 


20050612 


< 141 




20050711 


< 189 




20050903 


< 198 


424 


20050612 


< 225 




20050709 


< 180 




20050711 


< 192 




20050903 


< 204 


436 


20050612 


< 138 




20050722 


< 156 


445 


20050612 


< 171 




20050903 


< 225 


450 


20050612 


< 132 




20050903 


< 189 


451B 


20050816 


< 393 




20050824 


< 243 


473AB 


20050824 


< 180 


486 


20050612 


< 129 


493.1 


20050824 


< 183 


504 


20050612 


< 159 




20050722 


< 144 


514 


20050612 


< 117 




20050722 


< 132 


526 


20050612 


< 177 




20050722 


< 153 


566A 


20050609 


< 141 




20050722 


< 144 
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Table 3 — Continued 



GJ Epoch Flux Density Limit 

(YYYYMMDD) (/xJy) 



569AB 


20050609 


< 147 


570B 


20050824 


< 285 




20050824 


< 285 


581 


20050612 


< 117 


625 


20050612 


< 129 




20050711 


< 153 


628 


20050816 


< 252 




20050824 


< 192 


638 


20050612 


< 120 




20050711 


< 153 


643 


20050816 


< 189 




20050824 


< 171 


644AB 


20050816 


< 171 




20050824 


< 171 


673 


20050609 


< 111 




20050711 


< 156 


687 


20050612 


< 126 




20050709 


< 162 


694 


20050612 


< 129 




20050711 


< 150 


695BC 


20050816 


< 270 




20050822 


< 240 


699 


20050609 


< 123 




20050711 


< 153 


701 


20050609 


< 123 




20050711 


< 171 


725A 


20050612 


< 141 




20050709 


< 168 


725B 


20050612 


< 141 




20050709 


< 168 
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Table 3 — Continued 



GJ 


Epoch 


Flux Density Limit 




(YYYYMMDD) 


/ T \ 

(/^Jy) 


747R 


90050809 


<r AA^ 

\ -L -L X 










90050609 


^ 144 




2005071 1 


< 171 

X 1 X 


7Q1 9 


90050899 


^910 

\ Zj XU 




90050fi1 9 


^ 19'^ 




900^070Q 




OUc/ 


90050fi19 






90050709 


^ 189 




90050609 


<^ 144 




9005071 1 


< 967 


A 


90050895 


< 931 






<r 97 




90050fi1 9 


< 97*^ 




9005061 2 


< 916 




2005070Q 


< 995 


OU / -rt. 




<■ 109 


87fiA 


20050609 


<^ 919 




20050612 


<^ 901 




20050709 


<^ 186 


ooU 


zuuououy 


^ 1 on 




20050612 


< 114 




20050709 


< 162 


896B 


20050807 


< 162 




20050825 


< 186 


905 


20050824 


< 174 


908 


20050612 


< 105 




20050709 


< 123 


1156 


20050824 


< 243 


1245ABC 


20050809 


< 219 
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Table 3 — Continued 



GJ 


Epoch 


Flux Density Limit 




(YYYYMMDD) 


(/^Jy) 




20050825 


< 192 


1289 


20050824 


< 459 


2005BC 


20050807 


< 1254 




20050825 


< 420 


3076 


20050807 


< 537 




20050825 


< 264 


3125 


20050816 


< 414 


3193B 


20050809 


< 198 


3304 


20050809 


< 339 




20050824 


< 249 


3820 


20050824 


< 207 


4053 


20050816 


< 159 




20050822 


< 150 


4274 


20050824 


< 198 
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Tablc 4. VLBA Observations of M dwarfs 





jjate 


Beam Sizes 


Beam PA 


T-i^<-i rr/^ "D1\/TQ 

image rvivio 






(mas, mas) 


(deg) 


/iJy 


65B 


23 MAR 2006 A 


( 3.5, 1.5) 


-10.7 


163 


65B 


23 MAR 2006 B 


( 4.2, 1.8) 


11.4 


141 


65B 


25 MAR 2006 A 


( 3.4, 1.4) 


-15.7 


157 


65B 


25 MAR 2006 B 


( 2.8, 2.0) 


8.4 


118 


65B 


26 MAR 2006 A 


( 3.5, 1.5) 


-13.6 


147 


65B 


26 MAR 2006 B 


( 2.7, 2.1) 


9.1 


121 


285 


21 MAR 2006 


( 2.1, 0.9) 


-1.1 


90 


285 


24 MAR 2006 


( 2.5, 1.0) 


13.4 


109 


285 


27 MAR 2006 


( 2.7, 1.1) 


15.9 


111 


412B 


25 MAR 2006 


( 2.4, 0.9) 


-11.2 


109 


412B 


30 MAR 2006 


( 2.2, 1.2) 


-12.7 


86 


412B 


01 APR 2006 


( 2.2, 0.9) 


-12.8 


86 


803 


10 MAY 2006 


( 3.2, 1.9) 


5.3 


111 


803 


20 MAY 2006 


( 3.2, 1.6) 


6.0 


111 


803 


21 MAY 2006 


( 3.3, 1.7) 


8.2 


109 


896A 


23 MAR 2006 


( 2.6, 1.5) 


2.3 


89 


896A 


25 MAR 2006 


( 2.4, 1.0) 


-10.9 


85 


896A 


26 MAR 2006 


( 2.4, 1.0) 


-10.2 


88 


1224 


10 MAY 2006 


( 2.9, 2.0) 


12.3 


88 


1224 


20 MAY 2006 


( 3.0, 1.6) 


-0.5 


104 


1224 


21 MAY 2006 


( 2.9, 1.8) 


0.6 


102 


4247 


23 MAR 2006 


( 2.4, 1.6) 


-5.2 


95 


4247 


25 MAR 2006 


( 2.1, 1.1) 


-16.6 


85 


4247 


26 MAR 2006 


( 2.2, 1.1) 


-16.2 


87 



Table 5. Positions of Astrometric Calibrators 



GJ 


Calibrator 


Flux 
(mJy) 












a 








6 


Type 


65B 


J0132-1654 


1294.0 ± 


5 


.7 


01 


32 


43. 


.487468 








-16 54 48.522043 


primary 




J0142-1633 


38.4 ± 


1 


.2 


01 


42 


45. 


,261198 


± 


00. 


.000014 


-16 33 07.196786 ± 00.000121 


secondary 


285 


J0739+0137 


463.6 ± 


2 


.2 


07 


39 


18. 


,033892 








01 37 04.617926 


primary 




J0751+0152 


50.4 ± 


1 


.4 


07 


51 


02. 


,281590 


± 


00. 


.000002 


01 52 15.761011 ± 00.000023 


secondary 


412B 


J1108+4330 


368.6 ± 


1 


.6 


11 


08 


23. 


,476932 








43 30 53.657076 


primary 




JlllO+4403 


103.7 ± 


2, 


.6 


11 


10 


46. 


,345812 


± 


00. 


,000003 


44 03 25.925528 ± 00.000028 


seconda^ 


803 


J2056-3208 


594.5 ± 


4 


.4 


20 


56 


25. 


.070234 








-32 08 47.800507 


primary 




J2042-3152 


62.6 ± 





.6 


20 


42 


46. 


.406548 


± 


00. 


.000075 


-31 52 28.592351 ± 00.002223 


secondary 


896A 


J2328+1929 


119.4 ± 





.8 


23 


28 


24. 


.874755 








19 29 58.030041 


primary 




J2334+2010 


62.2 ± 


1 


.9 


23 


34 


14. 


.156496 


± 


00. 


.000009 


20 10 28.882640 ± 00.000279 


secondary 


1224 


J1753-1843 


110.4 ± 





.5 


17 


53 


09. 


,088754 








-18 43 38.523184 


primary 




J1825-1718 


69.1 ± 





.3 


18 


25 


36. 


,532398 


± 


00. 


,000029 


-17 18 49.849551 ± 00.002606 


secondary 




J1809-1520 


16.4 ± 





.2 


18 


09 


10. 


,208891 


± 


00. 


,001145 


-15 20 09.692410 ± 00.030386 


secondary 


4247 


J2205+2926 


141.5 ± 





.6 


22 


05 


46. 


,506426 








29 26 55.131163 


primary 




J2203+2811 


56.8 ± 


1 


.5 


22 


03 


59. 


,147970 


± 


00. 


,000004 


28 11 21.869208 ± 00.000015 


secondary 



Table 6. Positions of Stars 



GJ 


MJD 


Flux (mjy) 


RA 


Dec. 


Aa 


AS 


Aaopt 


A5opt 








(J2000) 


(J2000) 


(mas) 


(mas) 


(mas) 


(mas) 



65B 



53817.78 
53817.89 
53819.89 

53820.89 



2.171 ± 0.482 
1.800 ± 0.340 
1.872 ± 0.289 

2.075 ± 0.443 



01 39 02.990153 ± 0.000014 
01 39 02.990308 ± 0.000011 
01 39 02.992267 ± 0.000008 

01 39 02.993345 ± 0.000016 



-17 56 57.918175 ± 00.000342 
-17 56 57.918094 ± 00.000287 
-17 56 57.906928 ± 00.000202 

-17 56 57.902321 ± 00.000456 



0.000 ± 0.198 
2.208 ± 0.161 
30.158 ± 0.116 

45.551 ± 0.229 



0.000 ± 0.342 
0.081 ± 0.287 
11.247 ± 0.202 

15.854 ± 0.456 



319.940 ± 155.000 
320.487 ± 155.000 
319.378 ± 155.000 

320.182 ± 155.000 



264.745 ± 155.000 
264.289 ± 155.000 
266.167 ± 155.000 

266.183 ±1155.001 
'~i;35.409 
ij: 35.409 
± 19.137 
± 28.187 
± 28.187 
± 28.186 
"85.510 ± 90.706 
85.268 ± 90.706 
85.804 ± 90.706 



285 53815.14 
53818.14 



1.476 ± 0.275 
0.553 ± 0.100 



07 44 40.017980 ± 0.000006 
07 44 40.017550 ± 0.000007 



03 33 6.109089 ± 00.000163 
03 33 6.108662 ± 00.000208 



0.000 ± 0.088 
-6.436 ± 0.112 



0.000 ± 0.163 
-0.427 ± 0.208 



-45.984 ± 41.247 
-46.170 ± 41.247 



73.853 I 
_74.115 
~-3.153 
"~20.565 
20.924 
20.572 



803 53865.20 1.772 ± 0.310 20 45 09.679701 ± 0.000011 -31 20 29.507331 ± 00.000290 



0.000 ± 0.135 



0.000 ± 0.290 



17.516 ± 26.115 



896A 



53817.81 

53819.81 
53820.81 



0.929 ± 0.203 

1.399 ± 0.298 
3.895 ± 0.230 



23 31 52.433685 ± 0.000008 
23 31 52.434263 ± 0.000013 
23 31 52.434540 ± 0.000002 



19 56 13.712351 ± 00.000171 

19 56 13.714736 ± 00.000202 
19 56 13.715419 ± 00.000045 



0.000 ± 0.112 

8.148 ± 0.183 
12.059 ± 0.031 



0.000 ± 0.171 

2.385 ± 0.202 
3.068 ± 0.045 



114.239 ± 35.428 

114.366 ± 35.428 
114.281 ± 35.428 



4247 



53817.79 
53819.79 
53820.79 



0.950 ± 0.201 
1.161 ± 0.300 
1.164 ± 0.257 



22 01 13.304385 ± 0.000008 
22 01 13.304782 ± 0.000016 
22 01 13.304933 ± 0.000011 



28 18 25.130036 ± 00.000173 
28 18 25.132249 ± 00.000257 
28 18 25.134030 ± 00.000196 



0.000 ± 0.104 
5.247 ± 0.207 
7.242 ± 0.148 



0.000 ± 0.173 
2.213 ± 0.257 
3.994 ± 0.196 



7.147 ± 122.891 
7.392 it 122.892 
6.914 ± 122.891 
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Table 7. Proper Motion and Acceleration Relative to Optical Astrometry 



GJ 








as 




(masy~^) 


(masy~^) 


(AU y-2) 


(AU y-2) 


65B 


37.3 ± 52.5 


-218.5 ± 38.8 


0.0032 ± 0.0045 


-0.0189 ± 0.0033 


285 


22.4 ± 12.5 


-32.1 ± 22.7 


0.0088 ± 0.0049 


-0.0127 ± 0.0090 


896A 


-7.4 ± 7.2 


-10.1 ± 24.5 


-0.0031 ± 0.0030 


-0.0042 ± 0.0102 


4247 


18.1 ± 27.8 


-24.4 ± 30.3 


0.0046 ± 0.0071 


-0.0062 ± 0.0078 



Table 8. Limits in Companion Mass and Semi-Major Axis 



GJ 

' min ,min 

(AU) (Mj) (Mj) (Mj) 



65B 


1.3 


4.5 


19 


24 


285 


0.8 


3.3 


9 


42 


896A 


0.8 


3.8 


10 


42 


1217 


1.0 


6.1 


21 


11 



